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Recently the CRESST collaboration has published the long anticipated results of their direct 
Dark Matter (DM) detection experiment with a CaW04 target. The number of observed events 
exceeds known backgrounds at more than 4(t significance, and this excess could potentially be due 
r~i to DM scattering. We confront this interpretation with null results from other direct detection 

experiments for a number of theoretical models, and find that consistency is achieved in non- 
minimal models such as inelastic DM and isospin- violating DM. In both cases mild tension with 
^ constraints remain. The CRESST data can, however, not be reconciled with the null results and 

I— I with the positive signals from DAMA and CoGeNT simultaneously in any of the models we study. 
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1. INTRODUCTION 



At the recent TAUP 2011 conference the CRESST collaboration has reported a significant 
(> 4a) excess of nuclear recoil events in their detector, which cannot be explained by any 
known source of background [1]. If interpreted in terms of dark matter-nucleus scattering, 
these events would point to dark matter (DM) masses roughly between 10 GeV and 60 GeV, 
and DM-nucleon scattering cross section of order 10~^^-10~'^° cm^. The CRESST result thus 
provides a third inconclusive hint for the existence of relatively light dark matter, in addition 
to the long-standing observation of an annual modulation signal by DAMA [2], and the more 
recently reported event excess and annual modulation signal from CoGeNT [3]. 

It is, however, well known [1, 4-7] that the DAMA, CoGeNT and CRESST results, 
when interpreted in terms of the simplest dark matter models, are inconsistent among each 
other, and are moreover in severe confiict with null results from other experiments, such as 
CDMS [8, 9], XENON-100 [10] and KIMS [11, 12] (for recent phenomenological analyses 
see Refs. [13-39]). The updated KIMS result [12] is of particular interest because KIMS' 
target material CsI(Tl) contains iodine, which is also one of the target elements used in 
the DAMA's Nal(Tl) crystals. The null result from KIMS rules out any explanation of the 
DAMA annual sinusoidal modulation signal in terms of iodine recoils. Thus, if DAMA is 
actually seeing dark matter, the scattering would have to be dominantly off sodium. 

In this paper, we focus on the CRESST result and discuss possible interpretations of the 
observed events in terms of dark matter scattering. We confirm and quantify that within the 
conventional scenario of elastic spin-independent scattering the parameter region required 
by CRESST is disfavored by constraints from other experiments. Therefore we focus on 
modified particle physics scenarios, and we show that inelastic scattering or isospin- violating 
interactions may alleviate the tension between CRESST and the null results, though a joint 
explanation with DAMA and/or CoGeNT seems not to be possible. Let us note that the 
naturally most abundant isotopes in CRESST's CaW04 target have zero spin, and therefore 
CRESST cannot efficiently probe spin-dependent dark matter interactions. 

The outline of the paper is as follows. In sec. 2 we review the relevant equations for 
DM direct detection (see also [4, 24, 40] for more details) and discuss the particle physics 
scenarios we are going to use to describe the data. In sec. 3 we list the data we are going 
to analyse and provide some details on how we do the fitting of the various experimental 
data sets. In sec. 4 we consider conventional elastic spin-independent scattering and confirm 
the severe tension between CRESST and other experiments (sec. 4.1). We also consider, in 
sec. 4.2, the implications of a (possibly large) contamination of the CoGeNT event excess due 
to inefficient rejection of background events occuring close to the surface of the detector [41]. 
In sec. 4.3 we discuss the possible impact of new preliminary measurements of the quenching 
factor for sodium recoils in DAMA's target material Nal(Tl) [41]. Next, we turn to more 
exotic particle physics models of DM, considering inelastic scattering (sec. 5) and isospin- 
violating interactions (sec. 6). In both cases we show that a description of CRESST data 
consistent with all null results is possible. We conclude in sec. 7. 
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2. DIRECT DETECTION AND DARK MATTER MODELS 



The differential rate for a DM particle x to scatter elastically in a detector composed of 
nuclei with charge Z and mass number A is given by 

dR 1 f 3 daA N /.X 

- ' dv—vf{^r). (1) 



dEd ruA J\v\>v^;^ dE^ 

The units in which we will quote dR/dE are events/keV/kg/day. Here is the local 
DM density for which we adopt the standard value of p^ = 0.3GeVcm~'^. All our cross 
sections are understood relative to this reference value. More realistic DM densities (see 
e.g. [42]) lead to a trivial re-scaling of our results, /(v) is the DM velocity distribution in 
the rest frame of the detector. It is obtained from the velocity distribution in the galactic 
rest frame, for which we assume a Maxwellian distribution with a smooth cut-off, /gai(v) oc 
[exp(— v^/t;^) — exp(fgg^/'{;^)]9(fesc — v), where v = 220kms~^ and Vesc = 550 kms""*^ (for 
more details see [24]). The lower limit of the integration in (1) is set by the minimal velocity 
fmin that the incoming DM particle has to have in order to be able to deposit an energy Ed 
in the detector. For the case of inelastic scattering, x'^i it is given by 



1 / itlaE, 
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+ 5 , (2) 



where 5 = m^i — m^ <^ m^. The same equation also applies to elastic scattering, with 5 = 0. 

In the scattering rate (1) the information about the particle physics DM model is encoded 
in the cross section daA/dEd- Since DM is nonrelativistic, the cross section can in general 
be divided into contributions from couplings to the mass of the nucleus — giving the spin 
independent (SI) cross section — and couplings to the spin of the nucleus — giving the spin 
dependent (SD) contribution to the cross section 

The SI cross section for DM-nucleus scattering can be written in terms of the SI cross section 
for DM-nnc/eon scattering: 

SI ^ [^/p + (^ - Z)fnf ^^lN SI (A\ 

f2 ,,2 V • ^ > 

Jp f 



XP 



Here, fp^n is the normalized SI DM coupling to protons and neutrons, respectively, p^^p is 
the reduced mass of the DM-nucleon system, and 0"^^ is the SI cross section for scattering 
of DM on a proton. For the SI form factor F{Ed) we use [43] F{Ed) = 3e~''^*^/^[sin(Kr) — 
Kr cos{Kr)] / (nr)^ , with s = 1 fm, r = \/R^ — 5s^, R = 1.2^4-'^/^ fm, k = y/2/mjqEd (and 

~ —H?). In this paper, we will not consider spin-dependent scattering (i.e. we set ct^° = 
in eq. (3)) since the isotopes which are present in the CRESST detector either have no 
spin, or have a very small natural abundance, so that CRESST is not very sensitive to 
spin-dependent DM interactions. 

We will analyze the following particle physics models for DM-nucleus interactions: 
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• Elastic spin independent (eSI) DM scattering. This type of DM-nucleus inter- 
action can for instance arise from excliange of a scalar particle (such as a Higgs boson 
or a squark), or from vector interactions (such as the vectorial part of Z exchange). 
It is thus very common in many models, examples include supersymmetric neutralino 
DM [43] , the lightest Kaluza-Klein particle in universal extra dimensions [44, 45] , and 
little Higgs models [46] . Because DM is nonrelativistic and the momentum transfer in 
DM-nucleus scattering is small compared to typical nuclear energy scales, eSI scatter- 
ing occurs coherently off all nucleons in a target nucleus, so that the cross section is 
enhanced by a factor (the effect of nonzero momentum exchange is accounted for 
by the nuclear form factors). Conventionally, eSI scattering is assumed to be isospin- 
conserving, fn = fp = 1- For instance, in models where DM scatters through Higgs 
exchange, this is because the process predominantly probes the sea quark content of 
the nucleon, which is the same for protons and neutrons. We will relax the assumption 
of isospin-independent scattering below. 

• Inelastic DM (iDM). In iDM models, the DM sector is composed of at least two 
particles with mass splitting 5 > 0, so that in a scattering process on a nucleus, 
the lighter state is converted into the heavier one. In this case, it is possible that DM 
scattering can give a signal above the typical (9(10 keV) detection thresholds only when 
the target nuclei are relatively heavy and the DM velocity is high, close to the galactic 
escape velocity fesc- This is required for the momentum exchange to be large enough to 
overcome the splitting [47]. iDM therefore avoids constraints from searches performed 
on lighter nuclei, and the annual modulation of the rate due to the change in relative 
velocity between the Earth and the galactic DM halo is enhanced. The possibility of 
exothermic scattering, 6 < 0, has also been discussed in the literature [36]. 

• Isospin-violating DM (IVDM) scattering. In general DM can couple differently 
to protons and neutrons, so that fp ^ fn [48]. Note that even DM interactions 
through Higgs exchange in supersymmetric models are isospin-violating because of 
the different up and down Yukawa couplings, and only appear isospin-conserving in 
direct detection experiments because couplings to sea quarks are dominant there. In 
the presence of new interactions, however, it is quite possible that isospin-violation in 
DM-nucleus interaction is much stronger. Note that isospin-violating DM does not 
imply large flavor-changing neutral currents. For instance, there could be different 
couplings of dark matter to left-handed and right-handed quark vector currents, see 
e.g. the interactions discussed in [49]. Isospin-violating DM can also be realized in 
technicolor models, see [50] . If DM couplings to protons and neutrons have opposite 
sign, fpfn < 0, there can be a cancellation between the two contributions, depending 
on the number of protons and neutrinos in a given target isotope. For example, 
for fn/ fp — —0.7 the rate of scattering on Xe is reduced by an order of magnitude 
[37, 48, 51]. We will investigate whether it is possible to use this effect to re-concile 
the CRESST result with the constraints from other experiments (see also [4, 5, 20] for 
recent analyses in the context of CoGeNT). 

If the DM interaction is mediated by a particle lighter than the typical momentum trans- 
fers in direct detection experiments, of order 100 MeV, the expected energy spectrum will 
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be modified. Recent studies [4, 5, 15] have shown that this does not improve the consistency 
of the DM hints with various constraints and therefore we will not discuss light mediators 
here. 

3. FITTING DIRECT DETECTION DATA 

In this section, we discuss the inputs to our fits for the various direct detection experiments 
we consider. 

3.1. CRESST 

CRESST-II, which uses CaW04 crystals as target material, obtains after 739 kg days of 
exposure 67 events that fulfill all acceptance criteria [1]. In order to reject backgrounds, 
the CRESST detectors record for each event both the energy deposited in phonons and the 
amount of scintillation light emitted. The scintillation light yield can be used to distinguish 
the signal from most sources of background, and also to distinguish scattering processes on 
the different target nuclei, W, Ca and O. The known backgrounds in CRESST are electrons 
and gamma rays from radioactive decays, a particles, neutrons, and ^^°Pb decays due to 
surface contamination. The total deposited energy in each event is measured in the phonon 
channel, where the lower energy threshold is chosen separately for each of the eight detector 
modules (see Table 1 in [1]), according to the requirement that the expected number of 
electro-magnetic background events in the signal window is less than one. The CRESST 
collaboration has found that only about 65-75% of the observed events can be accounted 
for by the known backgrounds, corresponding to an excess with a statistical significance of 
more than 4ct. 

While the CRESST collaboration performs a likelihood fit of DM parameters to their 
data, using information on the energy deposit and the light yield for each event, we perform 
a simplified fit utilizing only publicly available information, which does not include the light 
yield for each event. In particular, we fit the total event rate in each detector module, as 
well as the overall energy spectrum. We take into account the individual energy threshold 
of each detector module. Following the CRESST collaboration, we assume backgrounds due 
to Pb recoils, a particles and neutrons to be constant in energy (see fig. 11 in [1]), whereas 
the e/7 background, which is one event per detector module, is assumed to appear in the 
lowest energy bin of each module. We compare our predicted energy spectrum to the data 
shown in fig. 11 of [1] using a analysis, in which the normalization of the flat background 
is allowed to float independently for each detector module. 

3.2. CRESST commissioning run 

In an earlier commissioning run of the CRESST experiment [52] using only two detector 
modules, the collaboration has collected a combined exposure of 47.9 kg days and observed 
three candidate events. The energy threshold of 10 keV in that analysis was lower than in 
the full run, whereas the upper end of the signal windows at 40 keV was the same. Using 
updated acceptance factors and considering DM scattering on all three target elements (the 



5 



CRESST collaboration only considered tungsten recoils in that analysis), tight constraints 
are obtained, excluding part of the preferred region from the longer CRESST-II run [53]. 
We include this constraint in our analysis using the energy and isotope dependent efficiences 
given in [53], and we apply the maximum-gap method [54] to compute an exclusion limit. 

3.3. CoGeNT 

The CoGeNT collaboration has first reported an excess of low-energy events — possibly 
due to dark matter — in 2010 [55]. Subsequently, this excess has been confirmed, and it 
has been found to show a seasonal variation at 2.8a significance [3]. Very recently, new 
investigations of CoGeNT 's backgrounds have indicated the possibility that a yet unknown 
fraction of the observed excess events is due to activity close to the surface of the detector, 
where external radioactivity is prevalent. Our analysis of the CoGeNT data from [55, 56] 
is similar to the one presented in [6]. We carry out an unbinned extended maximum likeli- 
hood fit, thus making the best possible use of the recorded energy of each individual event. 
In addition to the peaked cosmogenic backgrounds and the flat background considered in 
ref. [6], we also include an exponential component of the form aex]i{—E/Eo) to account for 
the possible contamination of the event sample by surface events. Since the values of the 
parameters a and Eq are not known precisely yet, we will consider different possibilities. 
We do not include timing information in our fit, since it is conceivable that not only the 
dark matter signal, but also the surface backgrounds show seasonal variation — for instance 
due to fluctuating levels of environmental radioactivity — and it would be very difficult to 
disentangle these different sources of annual modulation, especially with the still relatively 
low statistics in CoGeNT. 

3.4. XENON-100 

In a 4843 kg day exposure of the XENON-100 detector, three candidate events have been 
observed, with a background expectation of 1.8 ± 0.6 events [10], so there is no evidence for 
DM scattering. These data lead not only to the world's strongest limit on conventional spin- 
independent elastic WIMP-nucleon scattering, constraining the corresponding cross section 
to below 10~^^ cm^ for ~ 100 GeV, but have also important implications for the low 
mass region around 10 GeV (see also the discussion in [57-65]). We derive a conservative 
limit by using the maximum gap method from [54] . 

An important input to the limits from xenon experiments is the scintillation light yield, 
parametrized conventionally by the function Lg^^E^^), see e.g. [34]. We use for Leff(-Enr) 
the black solid line from fig. 1 of [10]. It has been obtained from a fit to various data, 
but is dominated at low energies by the recent measurement from ref. [66]. We include an 
overall systematic uncertainty on L^s, which we assume to be ±0.01 above E^^ = 3 keV, and 
increasing linearly by 0.05/keV below. The assumptions on Leff below Ej^^ ^ 3 keV (where 
no data are availble) have no impact on the XENON-100 exclusion curves. Even assuming 
Leff = below 3 keV, the limits remain essentially unchanged [4]. Hence, the result does 
not rely on any extrapolation into a region with no data. 
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3.5. KIMS 



The KIMS experiment uses a Csl target, and is thus directly probing whether the DAMA 
signal can be due to DM scattering on iodine. The KIMS collaboration have recently pre- 
sented results based on 24524.3 kg days of exposure, corresponding to roughly 1 year of 
data taking [12]. Since they do not see any events at the relevant recoil energies below 
8 keVee^, they are able to exclude the possibility that the DAMA annual modulation can 
be due to WIMPs recoiling on iodine. Note that this statement is very robust with respect 
to astrophysical uncertainties and our ignorance of the underlying particle physics model of 
DM. KIMS results also place strong constraints on inelastic DM because both Cs and I are 
heavy, so that iDM scattering is very efficient in KIMS. 

3.6. CDMS-II 

CDMS-II bounds are derived from the 612 kg days of data taken with Ge detectors 
in four periods between July 2007 and September 2008 [8], and from the previous CDMS 
search which had an exposure of 397.8 kg days, obtained between October 2006 and July 

2007 [67]. In [8] two events with recoil energies of 12.3 keVnr and 15.5 keVnr were observed 
in the 10-100 keVnr signal window. The expected backgrounds from misidentified surface 
events, cosmogenic backgrounds and neutrons were 0.8 ± 0.1 ± 0.2, 0.041q;q3 and 0.03-0.06, 
respectively. In our fits we use again Yellin's maximum gap method [54] to set limits. We 
take into account a linear efficiency drop from 32% at 20 keVnr to 25% at both 10 keVnr 
and 100 keVnr, and we use a constant energy resolution of 0.2 keV. 

3.7. CDMS-Si and CDMS low threshold analysis 

A stringent constraint on low mass DM is obtained from a modified reanalysis of the 
CDMS-II data collected in eight germanium detectors between October 2006 and September 

2008 [9]. In this analysis, the CDMS collaboration lowered the analysis energy threshold 
to 2 keVnr at the expense of allowing for a larger background contamination, and obtained 
and enhanced sensitivity to DM-nucleus scattering for WIMP masses below ~ 10 GeV. 
In our analysis of the low-threshold data, we use the Ax^ method, but include only those 
bins in which the predicted number of events is larger than the observed one. The CDMS 
collaboration has also collected 12 kg days of data on silicon targets [68]. No signal events 
above expected background were observed. Since silicon is lighter than germanium this 
places competitive bounds on light DM, especially in the case of isospin- violating DM models. 



The units keVee are used for "electron equivalent energy" , i.e. the amount of energy deposited in electronic 
excitations. Usually, this is only a fraction (given by the quenching factor) of the initial nuclear recoil 
energy, for which we use the unit kcVnr. 
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3.8. DAMA/LIBRA 



The combined exposure of the DAMA/LIBRA experiment and the previous DAMA/Nal 
experiment (DAMA for short) is 1.17 ton yr and spans 13 annual cycles [2]. The data show 
an annual modulation signal at 8.9a significance, with a peak around the end of May. This is 
consistent with the expectation from WIMP scattering, which predicts a peak around June 
2nd for standard assumptions on the DM halo. In the fit we use the signal region from 2 keVee 
to 8 keVee. Above this range the data are consistent with no modulation. Furthermore, no 
signal is predicted above this range from the DM models considered, so higher energy recoils 
give no additional constraint on the fit and can be ignored. Our analysis of the DAMA data 
is analogous to the one presented in [24, 26, 40]. Unless otherwise noted, we use quenching 
factors gNa = 0.3 ± 0.03 and qi = 0.09 ± 0.009 [69] and ignore channeling effects, which are 
expected to be small [70]. Note that we include independent 10% systematic uncertainties 
on gNa and gj. 

4. ELASTIC SPIN INDEPENDENT SCATTERING 
4.1. Fitting CRESST results 

We first assume that the signal in CRESST is due to DM scattering elastically with spin 
independent couplings (eSI). The resulting allowed regions in DM mass and DM-proton 
scattering cross section ap assuming equal couplings to protons and neutrons (/p = /„) are 
shown in fig. 1. We obtain our best fit point at = 12.5 GeV and ap = 2.7 x 10^^^ cm^, 
while the degenerate solution with somewhat higher mass and smaller cross section, = 29 
GeV and ap = 10~^^ cm^, gives a description of data that is only slightly worse by Ax^ = 0.3. 
The agreement between our regions and the ones obtained by the CRESST collaboration 
validates our simplified procedure. As in [1] we find two local minima of at masses 
around 12 GeV and 25 GeV. Since we cannot use the information on the scintillation light 
distribution our fit is slightly less constraining than the CRESST collaboration's, and the 
parameter regions corresponding to the two local minima are actually merged into one 
region at 90% confidence level. Note that for small DM masses, scattering on W does not 
contribute since the attainable recoil energies on such a heavy target are below threshold. 
For higher DM masses, W becomes more important; for instance, the fraction of tungsten 
recoils reaches 90% for m-^ ~ 50 GeV [1]. The CRESST region for low DM masses has some 
overlap with the region where scattering on Na can explain the modulation signal in DAMA, 
under the assumption of a quenching factor gNa = 0.3 ± 0.03 in DAMA. 

It has been noted in [1] and stressed in [53] that data from the CRESST commissioning 
run [52] imposes relevant constraints in the parameter region preferred by currect CRESST 
data. The limit from the commissioning run is shown in fig. 1 as a black long-dashed curve, 
which excludes part of the region preferred by latest data, but a consistent region remains 
even at 90% CL. Since both runs use very similar detectors (up to minor differences in 
acceptances, backgrounds and energy thresholds), it is very difficult to alleviate the tension 
by invoking exotic particle physics. Since the commissioning run took place between the 
end of March and the end of July 2007 it is conceivable that non-standard DM halos with 
strong seasonal variation in the scattering rate could reduce the tension. Note that the 
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FIG. 1: Constraints on elastic, spin-independent, isospin-conserving DM-nucleon scattering. We show the 
parameter regions preferred by the CRESST-II and DAMA signals (for CoGeNT see fig. 2), together with 
constraints from KIMS, CDMS (high threshold and low threshold analyses), XENON-100 and the CRESST 
commissioning run. 

recent CRESST analysis and the commissioning run data are based on different acceptance 
cuts, and a direct comparison might be subject to systematic uncertainties. 

In order to quantify agreement or disagreement between data sets, we use the parameter 
goodness of fit (PG) test [71]. This test is based on the function 

Xpg = Axi + Ax2 with Ax- =X?(globalbf)-x-,min, (5) 

where the index i = 1,2 labels the data sets whose compatibility is to be tested, and 
is the difference between the of the i-th data set at the global best fit point (i.e., 
at the minimum of Xi + X2) the minimum x^ from a fit to the i-th data set alone. 
XpG measures the "price" one has to pay for combining the data sets, compared to fitting 
them independently. Xpg follows a x^ distribution with a number of degrees of freedom 
corresponding to the number of parameters to which both data sets are sensitive (see [71] 
for a precise definition). As shown in Table I, the PG test finds consistency between the 
full CRESST-II data set and the data from the commissioning run at the level of 10%. The 
combined best fit point is obtained at = 12.9 GeV and cxp = 2.0 x 10"^^^ cm^. 

For comparison we show in fig. 1 also constraints imposed on the eSI DM mass and cross 
section by various null searches, confirming that an interpretation of CRESST data in terms 
of elastically scattering spin-independent and isospin-conserving dark matter is ruled out by 
XENON-100 [10], CDMS [8], and the CDMS low threshold analysis [9]. As we can see from 
Table I, the PG test gives a probability for consistency between CRESST versus CDMS and 
XENON of less than 10"^ 

Below,we discuss modified particle physics models with the aim of bringing CRESST 
results into agreement with those bounds. Before we do that, however, let us briefly address 
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Data set 


A. mm 


[GeV] 


a[cm2] 5 [keV] U/ U ^x' [xld 


PC 


Elastic spin-independent scattering (eSI) 


CRESST 


27.7 


12.5 


2.7 X 10"^^ 1 


0.3 




CRESS Tcoi]^ 








1 


4.2 




Combination 


32.2 


13.2 


2.0 X lO"''^ 1 


[4.5] 


10% 


CRESST 


27.7 


12.5 


2.7 X 10"^^ 1 


25.8 




XelOO+CDMSLT+CRESSTcom 





1 


0.00 




Combination 


53.5 


0.8 


7.9 X 10-^5 1 


[25.8] 


< 10-5 


Inelastic spin-independent scattering (iDM) 


CRESST+CRESSTcom 


32.2 


13.2 


2.0 X 10"^^ 1 


7.6 




XelOO+CDMSLT 








1 


1.9 




Combination 


41.9 


34.9 


1.7 X 10-38 94 1 


[9.5] 


2% 


Isospin-violating dark matter (IVDM) 


CRESST+CRESSTcom 


31.9 


14.2 


3.9 X 10--^° -0.66 


4.1 




XelOO+CDMS+CDMSsi 








1 


9.2 




Combination 


45.1 


14.2 


2.6 X 10"^° -0.70 


[13.3] 


0.4% 



TABLE I: Consistency of different experimental data sets for various particle physics models using the 
parameter goodness-of-fit test [71] . We list the minimal for the different data sets and for their combina- 
tions, as well as the dark matter masses and dark matter-proton scattering cross sections at these minima. 
The column labeled "Ax^ [XpgI" lists, for the indvidual data sets, the Ax^ of the combined best fit point 
with respect to the individual best fit points, and Xpg = ^Xi + ^xi foi" the combined analyses. In the col- 
umn "PG" we give the probability of consistency from the parameter goodness of fit test, which is obtained 
by evaluating the cumulative distribution function of the distribution with the appropriate number of 
degrees of freedom (2 for eSI and 3 for iDM and IVDM, see text for details) at Xpg- 

recent developements concerning CoGeNT and DAMA. 
4.2. Increased background in CoGeNT? 

A recent reanalysis of CoGeNT data has revealed the possibility that a significant fraction 
of the so far unexplained low energy event excess is due to contamination by background 
activity on the surface of the detector ("surface events") [41]. We investigate in fig. 2 (left) 
what the implications of this increased background are for the deduced properties of DM 
(assuming the remaining signal is due to elastic, isospin-conserving, spin-independent dark 
matter-nucleus scattering). We parameterize the surface background as aexp{—E/Eo), 
where Eq is kept fixed at a value of 0.3 keVee, as suggested by the plots in [41]. Since 
the rate of surface events in the signal sample has not been determined precisely yet, we 
show results for a = keV"^ day~^ (red), a = 12 keV~^ day~^ (blue), and a = 24 keV~^ 
day~^ (green). At the respective best fit points in and dp, these values correspond 
to signal fractions in the 0.5-1.0 keVee interval of about 73%, 46% and 30%, respectively. 
Fig. 2 (left) shows that, as expected, larger background contamination shifts the CoGeNT- 
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WIMP mass m ^ [GeV] WIMP mass m ^ [GeV] 



FIG. 2: Left: Variation in the CoGeNT- favored parameter region for elastic spin-independent dark matter- 
nucleon scattering for different assumptions on the surface background in CoGeNT (see text for details). 
Right: The DAMA preferred regions for the standard assumption on the sodium quenching factor, q-^g_ = 
0.3, (orange region), for the smaller energy dependent gNa obtained in [41] (light red region), and for the 
hypothetical case where scattering on iodine is forbidden (dark red region). In all cases, we have assumed 
a 10% overall systematic uncertainty on the DAMA quenching factors. 

preferred region in the m^-a plane towards lower cross sections and increases its overall size. 
Similar conclusions have been reached in [41]. 

4.3. New sodium quenching factors for DAMA? 

The quenching factors in the DAMA experiment, which are required to convert the ob- 
served energy deposit (in units of keVee) to the true nuclear recoil energy from a dark matter 
interaction (in keVnr), have been under active discussion for some time already [70, 72, 73], 
and a new facet has been added to this discussion recently, when Collar et al. carried out 
an independent measurement of gNa, the quenching factor for sodium recoils in Nal(Tl) [41]. 
Their results, which are in tension with results from other groups [74-78], indicate that gNa 
may be lower than the standard value gNa = 0.3 used by the DAMA collaboration and in 
many phenomenological analysis. In particular. Collar et al. find gNa — 0.1 at 30 keVnr 
nuclear recoil energy and gNa — 0.2 at 200 keVnr. This pronounced energy dependence of 
gNa has also not been seen in previous measurements. 

In fig. 2 (right), we explore the implications that Collar et al.'s results, should they 
be confirmed, would have on the parameter space of elastic spin-independent DM-nucleon 
scattering. The orange regions show the parameter values preferred by DAMA at the 90% 
and 3a confidence levels for the standard assumption gNa = 0.3 ± 0.03, gi = 0.09 ± 0.009. 
(As before, we include a 10% systematic uncertainty on the quenching factors.) In this 
parameter region, the signal in DAMA is dominated by Na recoils since only a small fraction 
of (9(10 GeV) dark matter particles can transfer enough energy to an iodine nucleus to induce 
an event above the experimental energy threshold. The light red contours in fig. 2 (right) 
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show the preferred region one obtains when the sodium quenching factors are taken from [41]. 
Note that at the relevant DM masses ~ 60 GeV, the signal in DAMA would be dominated by 
iodine recoils. For comparison, we also show, in dark red, where the DAMA-allowed region 
would lie in a hypothetical model where iodine recoils are switched off. (Such a scenario 
could be realized, for instance, in the framework of resonant DM [79] or isospin-violating 
DM [51].) We see that in this case, the cross section required to explain DAMA remains 
roughly the same as with the standard quenching factors, but the preferred dark matter 
masses are larger, since more energy is required for an event to be above threshold. We 
conclude that, if the low values for g^a reported in [41] were confirmed, an interpretation 
of the DAMA result in terms of eSI scattering would become even more difficult than it 
already is for gNa = 0.3, since the corresponding allowed region is strongly excluded by 
contraints from XENON-100, CDMS, KIMS, and others. Note in particular, that there is 
no longer any region where the DAMA modulation is explained in terms of Na recoils. Iodine 
recoils, on the other hand, are disfavored in a rather model-independent way by KIMS [12]. 
For the DAMA regions shown in the remaining figures of this work we stick to the value 
gNa = 0.3 ± 0.03. 

5. INELASTIC SCATTERING 

We next investigate the possibility that the CRESST-II signal is due to inelastic scattering 
of DM on nuclei. Fig. 3 (left) shows the of CRESST as well as CRESST combined with 
XENON and CDMS as a function of the mass splitting 6. We observe that CRESST alone 
prefers elastic scattering, although a reasonable fit is possible up to 5 ~ 95 keV. However, 
the combined shows a pronounced minimum around 5 ~ 94 keV. We show the allowed 
region for CRESST and the exclusion limits from XENON-100 and KIMS in fig. 3 (right) 
for (5 ~ 90 keV. We observe that a CRESST region consistent with all constraints appears. 
The consistency according to the PC test has a probability of about 2% (see table I), which 
is a huge improvement compared to the elastic case. The still somewhat small value of 
this probability indicates that some tension is left between CRESST and limits from other 
experiments. 

Inelastic scattering favours heavy target nuclei. Using the very heavy tungsten nucleus 
{A ^ 184) the rate in CRESST can be enhanced relative to the rate in all other experiments, 
which leads to the improved consistency. On the other hand, for increasing S the signal events 
in an iDM scenario are shifted to higher energies, and since the CRESST data are more 
concentrated at low energies, too high values of 6 are disfavored by the data, explaining the 
behaviour of the curves in fig. 3. In fig. 4 we show the predicted spectrum at a parameter 
point close to the iDM best fit point. We see that, although it provides a reasonable 
description of the data, iDM cannot very well account for all of the observed excess at low 
energy. The total number of predicted signal events at this parameter point is 12.2. In 
comparison we show also the predicted spectrum at the best fit point for elastic scattering 
(shaded histogram), where the number of predicted signal events is 24.8. 

In case of iDM with 5 ~ 90 keV, scattering in CRESST occurs exclusively on W since 
scattering on the light O and Ca nuclei is kinematically forbidden. At the best fit point 
for elastic spin-independent scattering, with = 12.5 GeV, on the other hand, we find 
a negligible contribution from W recoils and O (Ca) recoils accounting for 47% (53%) of 
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FIG. 3: Left: the total for CRESST (recent data from [1] combined with the constraint from the com- 
missioning run following [53]) and for CRESST plus data from XENON-100 [10] and CDMS low-threshold 
[9] as a function of the mass splitting S for inelastic DM. Right: iDM solution to CRESST for S = 90 GeV, 
together with constraints from the CRESST commissioning run, from KIMS, and from XENON-100. The 
DAMA preferred region for that S is also shown. 
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FIG. 4: The observed CRESST event spectrum compared to predictions for a few benchmark models. 
We show results for elastic SI scattering (shaded area), inelastic SI scattering (dashed blue histogram) 
and isospin- violating SI interactions (solid red histogram). The parameter values given in the legend are 
understood to be in GeV for m^, cm^ for the cross section, and keV for 5. Thin curves show the background 
only, while thick curves include both the signal and the background. 



the rate. While W recoil events from inelastic DM scattering do provide a reasonable fit to 
the Enr spectrum shown in fig. 4 the question remains whether also the scintillation light 
yield distribution is consistent with scattering exclusively on W. Unfortunately there is not 
enough information publicly available to check this explicitly. Let us note however that 
the favored parameter regions found by the CRESST collaboration in their fit of elastically 
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scattering DM [1] contain solutions that correspond to W-fractions as large as 90%. This 
suggests that because of the large overlap between the acceptance bands for W, O and Ca, 
scattering dominantly on W is consistent with the data, The final answer to this question 
can be obtained, however, only by a full fit to the 2- dimensional i^ji^-light yield distribution. 

6. ISOSPIN- VIOLATING SCATTERING 

Finally, we return to spin-independent elastic interactions, but relax the assumption of 
equal couplings of DM to protons and neutrons. That is, we will now assume fp ^ fn in 
eq. (4). Following [4] we define an effective nuclear mass number A^s such that the event 
rate from scattering on a given element is suppressed by Al^/A"^ relative to the case fp = fn' 



Here tan^ = fn/ fp and is the relative abundances of the i-th isotope in the target. Fig. 5 
(left) shows the suppression factor for various relevant elements. This plot shows that for 
fn/ fp = —0.7 it is possible to suppress the rate on Xe and to a lesser extent also that on 
Ge relative to O and Ca. Indeed, we can see from fig. 5 (right) that the combined from 
CRESST, XENON, and CDMS has a clear minimum at that value of fn/ fp- From the left 
panel we expect that the scattering rate in CRESST will be dominated by Ca and O, since 
W has a neutron to proton ratio similar to Xe, leading to a suppression of the W recoil rate 
at similar values of fn/ fp- On the other hand. Si has a proton to neutron ratio very similar 
to Ca, which implies that the bound from CDMS Si remains important. 

Fig. 6 shows the CRESST allowed region compared to the various limits assuming 
fn/fp = —0.7 (this value is chosen in order to minimize the predicted rate in Xenon- 
100, thus relaxing the otherwise strongest constraint on dark matter scattering). We ob- 
serve that part of the CRESST 3cr region remains consistent with XENON and CDMS 
limits at 90% CL. The combined best fit point is obtained for = 14.2 GeV and 
Ceff = (<7p + o"n)/2 = 2.6 X lO"''^ cm^. As expected the event rate in CRESST is domi- 
nated by the light elements (39% scattering on O, 61% on Ca) with a negligible contribution 
from W. The corresponding spectrum is shown in fig. 4 as the red histogram. However, 
although there is clear improvement with respect to the isospin-conserving case, the consis- 
tency according to the PC test is still rather low, at 0.4%, see table I. When performing a 
combined analysis of CRESST+CDMS+XENON we find closed regions only at 2(j, whereas 
at 3(7 no closed region appears, indicating also considerable tension in the data. 

Let us mention also that we do not expect further improvement by combining isospin- 
violating interactions with inelastic scattering. For IVDM with fn/fp = —0.7 scattering 
in the preferred parameter region occurs dominantely on O and Ca, since W recoils are 
suppressed because of the proton-to-neutron ratio ~ 0.7 in W. Inelastic scattering, on the 
other hand, occurs exclusively on W, since scattering on the lighter targets is kinematically 
forbidden. Therefore, iDM and IVDM are not compatible with each other as an explanation 
for the CRESST signal. 




(6) 
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FIG. 5: Left: the suppression factor for isospin- violating interactions relative to the iso-spin conserving 
case for various elements. A^f^ is defined in eq. (6). Right: The total for CRESST (recent data from [1] 
combined with the constraint from the commissioning run following [53]) and for CRESST plus data from 
XENON-100 [10] and CDMS (Ge [8] + Si [68] data) as a function of the ratio of the neutron and proton 
couplings. 
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FIG. 6: CRESST preferred region in dark matter mass and average scattering cross section = 
{<jp + (T„)/2 for elastically scattering isospin- violating dark matter with fn/ fp = —0.7, and constraints from 
other experiments. 
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7. CONCLUSIONS 



In this paper, we have studied the compatibihty of the results from various direct dark 
matter detection experiments, including in particular new results from CRESST, which 
show a 4(7 excess of events above known backgrounds. We have found that none of 
the considered dark matter models — elastic and inelastic spin-independent scattering as 
well as isospin-violating dark matter — can simulatenously explain the positive signals from 
CRESST, DAMA and CoGeNT, while satisfying constraints from other experiments. The 
CRESST signal alone, however, is consistent with exclusions from other experiments in 
the case of inelastically scattering DM and DM with isospin-violating couplings to visible 
matter, though mild tension remains. The parameter region preferred by current CRESST 
data and the null result from an earlier commissioning run of the experiment are consistent 
at the 10% level, as estimated using the parameter goodness of fit method (for standard 
assumptions on the DM velocity profile). 

We have also briefly commented on a recent reassessment of the CoGeNT data, which 
indicates that a yet unknown fraction of the excess events in CoGeNT is due to surface 
backgrounds, and on a new measurement of the Na quenching factors in DAMA, which 
indicates that they might be much smaller than previously thought. We have demonstrated 
that a lower signal rate in CoGeNT could make their results more consistent with the 
null results from XENON-100 and CDMS. The modified quenching factors in DAMA, if 
confirmed, would on the other hand increase the tension between DAMA and the null 
results. 

We conclude that the landscape of dark matter direct detection is evolving rapidly, and 
there is currently no known way of explaining all the positive signals (DAMA, CoGeNT and 
CRESST) simultaneously with all the null results. It is therefore crucial that these signals, 
but also the null results, are scrutinized carefully, and experimental cross-checks are carried 
out to test them as model-independently as possible, for instance by using the same target 
material in several experiments. 

Note added: During completion of this work we learned that the inelastic DM solution 
for CRESST has been independently obtained by N. Weiner [80]. 
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